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Abstract 

In this paper, we will prove a very general result of stability for 
perturbations of linear integrable Hamiltonian systems, and we will 
construct an example of instability showing that both our result and 
our example are optimal. Moreover, in the same spirit as the notion of 
KAM stable integrable Hamiltonians, we will introduce a notion of ef- 
fectively stable integrable Hamiltonians, conjecture a characterization 
of these Hamiltonians and show that our result prove this conjecture 
in the linear case. 

1 Introduction and results 

1. Let n > 2 be an integer, T" = M"/Z" and B = Bn be an open ball in 
of radius -R > 1 with respect to the supremum norm. We shall consider 
a near-integrable Hamiltonian of the form 

{H[e,i) = h{i) + f{e,i) 

\\f\<e«l 

where {9,1) G x 5 are action-angle coordinates for the integrable part h 
and / is a small perturbation in some suitable topology defined by a norm 
I . I . For simplicity, we shall restrict ourself to the analytic case, that is we 
assume that h and / are bounded and real-analytic on D = T" x i?, so that 
they have holomorphic extensions to some neighbourhood 

V„{D) = {{6 J) € (C"/Z") X C" I \I{e)\ < a, d{I,B) < a} 

for some cr > 0, where I denotes the imaginary part and d is the distance 
associated to the supremum norm. Then the norm of the perturbation 
I/I = l/U is defined by 

l/l<T = |/lco(y,{D)) = sup |/(2;)|. 
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In the absence of perturbation, that is when e is zero, all solutions {0{t), I{t)) 
of the corresponding Hamiltonian system are quasi-periodic and the action 
variables I{t) are integrals of motion, but in general this is no longer the 
case after perturbation. Basically, assuming some non-degeneracy condition 
on the integrable part h and some regularity on h and /, KAM theory is 
concerned with the persistence of many quasi-periodic solutions ( |Kol54] ). 
while Nekhoroshev theory deals with the variation of the action components 
of all solutions r |Nek77p . 

In the realm of KAM theory, one can introduce a notion of "KAM-stable" 
integrable Hamiltonian h for which any sufficiently small perturbation pos- 
sesses a set of positive Lebesgue measure of quasi-periodic solutions closed 
to the unperturbed ones, and such that the measure m(e) of the complement 
of this set satisfy lim£_j.o m{e) = 0. These KAM stable Hamiltonians have 
been characterized: they are exactly Riissmann non-degenerate Hamiltoni- 
ans, that is functions h : B ^ M. such that Vh{B) is not contained in any 
hyperplane of M". We refer to the nice survey |Sev03| for precise results and 
references. 

It is tempting to define a notion of "Nekhoroshev-stable" Hamiltonians 
and to try to characterize them, but first one has to come up with a precise 
definition. It is easy to see from the perturbative character of the Hamil- 
tonian that for all solutions {9{t),I{t)) with initial condition {9o,I()), one 
has 



Without further hypothesis on h, this cannot be improved. Indeed, following 
|Nek79] and |Nie06j . if the restriction of h to some affine hyperplane, whose 
direction is generated by integer vectors, has a non-isolated critical point, 
then there exist 6 > and an arbitrarily small perturbation of size e such 
that 



A simple example of this type will be given below. This prompts us to 
introduce the following two definitions. 

Definition 1. An integrable Hamiltonian h : B —^M is rationally steep if its 
restriction to any affine hyperplane of the form Iq + A, with Iq G B and A a 
linear subspace o/M" generated by integer vectors, has only isolated critical 
points. 

Definition 2. An integrable Hamiltonian h : B ^ M is effectively stable 
if for any f : T"- X B ^ R with \f\ < e, all solutions {9{t),I{t)) of the 
Hamiltonian system H = h + f starting at {9q,Iq) satisfy 




sup \I{t)-Io\>6. 



0<t<e-i 



lim 

£->0 



sup \I{t) - /ol = 0. 

0<\t\<e-'^ J 
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Hence an integrable Hamiltonian is effectively stable if one can ensure 
that the time of stability T(e), that is the maximal time during which the 
variation V{e) of the action of all solutions satisfy lim£_>.o V{e) = 0, is at 
least 1/e. Then how large one can choose T(e) will depend on more specific 
properties of the integrable Hamiltonian. 

The following conjecture was essentially made in |BoulOc| . 

Conjecture 1. Effectively stable Hamiltonians are exactly rationally steep 
Hamiltonians. 

As we already explained, if a Hamiltonian is effectively stable, then it has 
to be rationally steep, so that only the converse statement of this conjecture 
is of interest. In this paper we shall prove this conjecture in the specific 
case where the integrable Hamiltonian is linear. This is an important case, 
not only because of its own interest, but also because it is a crucial step 
in showing the conjecture in its full generality (this point will be briefly 
discussed at the end of the paper). 

2. So from now on we shall restrict to linear integrable Hamiltonians. 
Given a vector uj G \ {0}, we let 1{I) = co.I be the linear Hamiltonian 
with frequency uj and we consider 



Without loss of generality, we may assume that the system has been rescaled 
so that |a;| = 1, therefore, reordering the components of u if necessary, we 
can write 

cj = (l,ai, . . . ,a„_i) = (l,a), |aj| < 1, i G {1, . . . , n - 1}. 

Such "degenerate" integrable systems are never KAM stable, but the prob- 
lem of effective stability is of different nature as we will show. 

First, one can easily see that / being non-rationally steep is equivalent 
to UJ being resonant, that is there exists A; S \ {0} such that k.uj = 0. In 
this case, the Hamiltonian 



H{ej) = i{i) + f{e,i) 

\f\<e«l. 



(*) 



H{ej) 



uj.i + f{e), f{e) 



esin{k.6) 



gives rise to a system which can be easily integrated: 




^0 + tuj [Z"] 

Jo + cos(A;.0o)- 



So any solution starting at {9q,Iq) with k.9o 



satisfy 



sup \I{t) - Iq 

0<i<£-i 



k\ > 1. 
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Hence non-rationally steep linear Hamiltonians are indeed not effectively 
stable, and our conjecture in this case states that linear Hamiltonians with 
a non-resonant frequency are effectively stable. However, the only case for 
which results are known is when the frequency satisfies a classical Diophan- 
tine condition, and there much stronger stability properties hold true in the 
sense that the time T(e) is at least exponentially large with respect to the 
inverse of the size of the perturbation. These results will be recalled below. 
In this paper, we shall prove a more general stability result and we will 
construct an example showing that it is the best one can obtain. 

3. Let us state precisely our results. We shall denote by | . 1^ the distance 
to the integer lattice Z, that is \x\z = min^g^ \x — p\ for x G M. The 
frequency vector u = (1, ai, . . . , a„_i) = (l,a) being non-resonant, the 
function ^ = given by 

^{K) = max{|A;.a|^^ | k £ I/"-^, < |A;| < i^} , K 

is well-defined. It is obviously strictly increasing on N*, hence we can extend 
it (keeping the same notation) as a strictly increasing continuous function 
defined on [l,+oo). Then let us also define two additional functions 

K{x)=x^{x), A(x)=A"^(x), X > 1, 

which are also strictly increasing and continuous. 

Theorem 1.1. Let H he as in with uj non-resonant. Then there exist 
positive constants Eq, c, ci and C2 depending only on n, R, a and uj such that 
if £ < Eq, all solutions {6{t),I{t)) of H with Iq € satisfy the estimates 

IHi) - Io\ < ci6, \t\ < 6e'^ exp (c2A(ce~^)) . 

for any (A {c£-^))'^ < ciS < R/2. 

Depending on the growth of the function A, the exponential factor in 
the time of stability below might not be very large but choosing 6 = with 
6 > arbitrarily small, thanks to the factor 6e~^ the time of stability is at 
least 1/e and this proves our conjecture in the linear case. 

Corollary 1.2. For linear integrable Hamiltonians, rationally steep Hamil- 
tonians are effectively stable. 

We shall give an elementary proof of Theorem 11.11 in the case n = 2, 
using only one rational approximation and a one-phase averaging, in the 
same spirit as in |Loc92] . This method of proof can be easily extended for 
any n > 2, but in general it fails to give the best result. 

Hence for any n > 2, we shall use more classical techniques, namely 
general resonant normal forms as in |Pos93j and [DG96j, and so our proof 
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will not be essentially new. In the special case where the frequency oj is 
Diophantine, that is when there exist 7 > and r > n — 1 such that for all 
k G \ {0}, \k.a\i > 7|A;|-^, then, in Theorem we can choose 

^(x) = ■y'^x'', A{x) = 7"-^3;^+^, A(x) = (jx)^ 
and our result gives 

\I{t) - Io\ < ci5, \t\ < 5e~^ exp (c2{c-fe'^)^^ . 
1 

for any (c~ 7" e) 1+^ < ci6 < R/2. In this case, this is exactly the result 
obtained in [DG96] . In |Pos93] . there is a similar but less flexible result, since 
there one can only choose 5 ~ 1. Let us also note that this Diophantine case 
was first considered in |Fas90] . and the result there was even more flexible 
since one has the freedom to choose 5 equals, up to a multiplicative constant, 
to e, but the proof is restricted to the Diophantine case and fails without 
this assumption. 

In any cases, our next theorem shows precisely that the time of stability 
achieved in Theorem 11.11 is the best possible. 

Theorem 1.3. For any non-resonant vector u G W^, there exists a sequence 
ifj)jm* of analytic functions on Va{D), with \ fj\a- = £j ^ when j — )• +00, 
such that the system Hj = 1 + fj has orbits which satisfy 

\I{t) - Io\ = \t\ejexp (^-2crA(ceJ^)) . 

for a constant c depending only on R, a and u. 

This result says that for a specific arbitrarily small perturbation and for 
some solutions of the perturbed system, the inequality of Theorem 11.11 are 
in fact equality (up to constants depending only on n,R,a and oj), so that 
Theorem 11.11 cannot be improved. Conversely, Theorem 11.11 implies that the 
above theorem is the best possible. 

The construction of the example of instability in Theorem 11.31 is elemen- 
tary, and it will follow very naturally from our proof of the stability result 
in the case n = 2. It is also inspired by an example given in |Sev03] . 

Therefore, as far as one is interested in exponential stability (that is, 
r(e) is exponentially large with respect to e~^), then a polynomial growth 
of the function ^, which is nothing but a Diophantine condition, is both 
sufficient and necessary. To obtain a non-trivial polynomial time of stability 
such as £~'^ , r > 1, then an exponential growth of the function ^ is sufficient 
and necessary. 

4. The plan of the paper is the following: the proof of the above theorems 
will be given in the next section, and further comments on the results are 
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given in the last section. Throughout the paper, in order to avoid cum- 
bersome expressions, we will replace positive constants depending only on 
n, i?, a and uj with a dot. More precisely, an assertion of the form "there ex- 
ists a constant c > depending on the above parameters such that u < cv" 
will be simply replaced with "u <• v" , when the context is clear. 

2 Proof of the results 

1. We shall start by giving an elementary proof of Theorem 11.11 in the 
special case n = 2. 

Given a vector v = {l,p/q) € M^, where p/q is a non-zero rational 
number in its lowest term, we let 1^ be the linear Hamiltonian with frequency 
V. The normal form result that we shall need, which is due to Lochak- 
Neishtadt ( |LN92j . see also |Pos99j ). is just a one-phase averaging. 

Lemma 2.1. Consider a Hamiltonian H = 1^ + f^j defined on V^(-D) with 
\fv\a <•£, and assume that for some K > 1, 

Kqe<l. (1) 
Then there exists an analytic symplectic transformation 

CD : V^/2iD) ^ V,{D) 
with 1$ — Id|o-/2 <■ q^ such that 

H> = Ho^ = l^ + g + f'^ 
with = 0, and the estimates 

bL/2<-e, \fv\<T/2<£e.-^ 

holds true. 

This is a so-called resonant normal form, that is the perturbation has 
been reduced, up to an exponentially small term to its resonant part g 
which satisfies {g^lv}- Now a simple observation is that if the denominator 
q is large, then so is the size of any integer vector which is orthogonal to 
V. In such a case, the resonant term g itself is exponentially small. Let us 
state this as a lemma which complements the result above. 

Lemma 2.2. Under the previous hypotheses, assume also that q > K . Then 
9 = 9 + 9, 9=1 9, \9\a/2<£, \9'\a/4<-£e~'^- 
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Proof. Expanding g as a Fourier series g{9,I) = XliteZ" gk{I)e^^^^'^ ■, we can 
write g = g + g' with 

g'{e,I)= 9k{I)e''^''-' . 

fceZ"\{o} 

Now the condition that {g, Z^,} = is easily seen to be equivalent to gk{I) = 
if k.v 7^ 0, that is 

g'{ej)= Yl 9k{I)e^'^''-'. 

k.v=0, k^O 

Now take k = (ki, k2) £ 7? \ {0} such that k.v = 0, that is qki + pk2 = 0. 
First one can see that neither ki nor k2 is zero, otherwise this would imply 
that either p or g is zero. Therefore q divides pk2, but as q and p are 
relatively prime, then q divides k2 and hence \k\ > \k2\ > q > K. This 
means that 

\9'l/2 < bUU/2, 9>KiO,I) = 9k{I)e^'^'-'. 



\k\>K 



It is now a classical estimate that 



lfi'>xU/4 <• bl<7/2e '^<ee 

with implicit constants depending only on a. □ 

Part of the arguments given above works only for n = 2. However, for 
any n > 2, some other simple arguments along these lines can be used (but 
as we already said, this does not give the best result in general). Now we 
can prove Theorem 11.11 for n = 2. 

Proof of Theorem \l.ll n = 2. Here lo = (1,q;) £ with \a\ < 1. The size 
of the perturbation e > being given, we choose K such that 

K^{K) = -e~^, K = A{-e-'^) 

for some implicit constant to be chosen below. First, assuming e is suffi- 
ciently small, that is e<l, we can ensure that K > 1 and we can apply 
Dirichlet's box principle to approximate a by a rational number with de- 
nominator q < ^{K): we obtain a vector v = {l,p/q) G M^, p G Z \ {0}, 
such that 

\uj -v\ <q-^^{K)-\ l<g<^(K). 
Moreover, from the definition of ^, 

*(g)"^ < \qa-p\ < 

so that q > K. Hence 

\uj-v\< K-^^'(K)-i=-e. 
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Now our Hamiltonian H can be written as 



H — I + f — Iv + fvi fv — I — Iv + f, 

with |/i,|o-<-e. Moreover, since q < ^{K), then qe<K~^, and choosing 
properly our imphcit constant in the definition of K, we can ensure that 
condition ([1]) is met and hence we can apply Lemma I2.lt there exists an 
analytic symplectic transformation 

$:V;/2(Z))^K(D) 
with |<I> - Id|^/2 <• <• such that 

H' = Ho<!> = l, + g + f^ 
with {g, l^} = 0, and the estimates 

\g\a/2<£, \K\a/2<^e~'^ 

holds true. Moreover, since q > K, Lemma 12.21 can also be applied and 
therefore 

H' = Ho^ = l,+g + g' + f:^ = k+g + f 
with g integrable and 

l/U/4 < b'L/4 + l/^L/4<-ee--^. 

Now let us write {9, 1) = ^{6', I'), and take any 6 such that <6 < R/2. 
Since ly + g is integrable, the mean value theorem gives 

\l'{t) - 1[)\<6, \t\<6e~^e^. 

Now coming back to the original coordinates and using |<I> — Id|^/2 <' 
this gives 

\I{t) - Io\<6, \t\<6e-^e^, 

which, recalling that K = A(-e~^), is our statement. This completes the 
proof. □ 



2. Now for any number of degrees of freedom, the strategy we explained 
above can be easily extended and a result of stability can be obtained. How- 
ever, it requires to compare simultaneous and linear Diophantine approxi- 
mation, and it seems to us that, unless the frequency is badly approximable, 
it cannot give the best result from a quantitative point of view. 

Therefore we shall use different arguments to prove Theorem 1 1.1 1 for any 
n > 2, following the "dual" approach. Given any sub-lattice A of Z", we 
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shall say that a function g defined on V^(-D) is A-resonant if its Fourier 
expansion is of the form 

fceA 

Then, given any K > I and any A > 0, we shall say that a vector w G M" is 
(A, i^)-non resonant modulo A if 

VA;GZ"\A, \k\<K, \k.w\>X. 

Finally, given such a vector w, let us write Iwil) = wJ. The normal form 
result that we shall use here is due to Delshams-Gutierrez ( |DG96) ). 

Lemma 2.3. Consider a Hamiltonian H = 1^ + f defined on Va{D), with 
l/lo- <• e and assume that 

K\-^e<l. (2) 
Then there exists an analytic symplectic transformation 

with |<D — Id|o-/2 <■ A^-'^e such that 

H' = Ho^ = l^+g + f 
where g is A-resonant, and with the estimates 

\9\a/2<£, \f'\a/2<£e~'^. 

The above normal form result strictly contains Lemma 12.11 (the latter 
corresponds to the case where A has rank n — 1, and we can choose A = 
for any K > 1), as a consequence its proof is more involved. A similar 
normal form result is contained in [Pos93] . but for our purpose it is weaker 
since the distance to the identity of the normalizing transformation $ is 
bigger (it is only of order KX~^e). 

The proof of Theorem 11.11 is now straightforward. 

Proof of Theorem \l.l[ Here a; = (1, qi, . . . , a„_i), and this frequency is by 
definition (^'(Er)~^, i^)-non-resonant modulo {0}. The size of the perturba- 
tion e > being given, we choose K such that 

K^{K) = -e~^, K = Ai-e^'^) 

with a suitable implicit constant so that the condition ([2]) is fulfilled. Then 
the requirement that K > 1 gives a threshold e <■ 1 . Therefore we can apply 
Lemma 12.31 there exists an analytic symplectic transformation 

$:V;/2(D)^K(D) 
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with \^ -ld\^/2<-'^{K)£ = - K-'^ such that 

H' = Ho^ = l+g + f' 

where g is {0}-resonant, that is g = g, and with the estimates 

l9U/2<-e, l/U/2 <-ee-^. 

Now, as before, let us write (0,1) = ^{9', I'), and take any 6 such that 
<• (5 < R/2. Since / + ^ is integrable, the mean value theorem gives 

\l'{t) - I^\<6, \t\<6£-^e^. 

Now coming back to the original coordinates and using |$ — Id 1 0-/2 <■ K~^, 
this gives 

-/o|<-'^, \t\<6e-^e^ 
which ends the proof, since K = A{-e^^). □ 

3. Now let show that our stability result is optimal by proving Theorem 1 1.3 1 

Proof of Theorem \1.3[ Given uj = (l,ai, . . . ,q„_i), let us pick one compo- 
nent, say ai, and let us denote by {pj/qj)jm the sequence of the convergents 
of tti. The vector u being non-resonant, ai is irrational and this means that 
the sequence (gj)jeN is strictly increasing. From the classical estimates 

{qj + qj+iy^ < \qjai - Pj\ < qj^^, j € N, 

we deduce that qj+i < "^{qj) < qj + qj+i < 2gj+i. Now the perturbation fj 
will be of the form 

f,i9J) = fj{I) + ff{0), (9,1) G X B. 
First, we choose fj{I) = Vj.I — lj.I, where Vj = {l,pj/qj, 02, • • • , ctn-i)- As 
|ai -PjQj^l < {QjQj+iy^ < 2('7j^'(gj))"\ 

if we set 

£j = c{qj^{qj))-\ qj = A{ceJ^) 

with a suitable constant c depending on ai, R and a, we obtain \ fj\o- < £i/2- 
Then, if we define kj = {pj, —qj,0, . . . , 0) € Z", we choose 

fj{6) = Sjfij cos{kj.0), fij = q~^ exp(-2crgj). 

For 6 G C" with \1{9)\ < a, we have 

\cos{kj.9)\ < exp(2cr|A;j I) = exp(2crgj) < qj/2exp{2aqj) 
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and therefore \fj\a < as well. Hence \fj\a < ^j, and ej — when 
j —7- +00. Now we can write the Hamiltonian 

Hj{e, I) = h{I) + fj{e, I) = Vj.I + ejq-^ exp{-2aqj) cos(%.0) 

and as kj.Vj = 0, the associated system is easily integrated: 

r e{t) = eo + tvj [z"] 

I I{t) = Iq — tkjqJ^Ej exp{—2aqj) cos{kj.9o). 

Choosing any solution with initial condition (Oq, Iq) satisfying kj.Oo = 0, 
cos{kj.6o) = 1 and using the fact that \kj\qj^ = 1, we obtain 

\I{t) — IqI = \t\ej exp{—2aqj). 

Recalling that qj = A(ceJ^) this gives 

\I{t) -Io\ = \t\ej exp (^-2aA(ceTi)) 
and this concludes the proof. □ 

3 Further comments 

Let us conclude this paper by several remarks. 

1. First, let us mention that a perturbation of a linear Hamiltonian system, 
as we have considered, also describes the dynamics in the neighbourhood of 
an invariant, linearly stable, isotropic and reducible torus carrying a quasi- 
periodic motion. Our Theorem 11.11 directly applies to the case where the 
torus is of full dimension (that is, when the torus is Lagrangian) and gives 
a non-trivial result of stability. In the case of intermediate dimensions, our 
result should apply also but there one has to work not only with angle-action 
coordinates but also with Cartesian coordinates. Anyway, in the case where 
the torus has dimension zero, that is for an elliptic fixed point, Theorem ll.il 
is useless since in this situation the perturbation has a peculiar form and 
classical Birkhoff normal form estimates apply and give a much better result 
of stability. Now Theorem 1 1 . 3 1 also applies in this setting, except for elliptic 
fixed point since in this case there is no proper angle-action coordinates to 
work with (and in fact, for n = 2, Diophantine and even Brjuno elliptic 
fixed points are always Lyapounov stable, see |Riis02] and |FK09j ). When 
the dimension of the torus is at least one, our result applies and in the 
Diophantine case, without any further assumptions, results of exponential 
stability are the best one can expect. Assuming further conditions on the 
Birkhoff invariants, it was shown in |MG95j (see also [BoulOaj for more 
general results) that super-exponential stability hold true. Here our result 
shows that some further conditions are indeed necessary. 
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2. Then, let us discuss the case where the system is assumed to be less 
regular than analytic. Our proof of Theorem 11.11 in the case n = 2 only 
uses a one-phase averaging. Such a result being available in the Gevrey 
category (|MS02j) and in the finitely differentiable category ( [BoulOb] ). the 
proof immediately extends in this case (of course, in finite differentiability, 
only polynomial stability hold true). Now for any n > 2, we have used a 
more elaborate normal form which is not known if the system is not ana- 
lytic. However, such a normal form can certainly be proved using smoothing 
techniques so that Theorem 11.11 should hold in any regularity. Obviously, 
the proof of Theorem 11.31 is independent of the regularity and the extension 
is straightforward. 

3. Finally, let us explain how our result can be used to prove in general that 
effectively stable Hamiltonians are exactly rationally steep Hamiltonians. As 
we have explained, our result deals with the case where the Hamiltonian is 
linear. In the non-linear case, first one has to use not only integrable normal 
forms but also more general resonant forms, which are also well-known. 
Then, the condition of being rationally steep roughly says that resonances 
"do not accumulate" so that any solution, when evolving (of course, if they 
do not evolve there is nothing to prove), will necessarily encounters non- 
resonant points, arbitrarily close and in any direction, and therefore the 
inductive scheme introduced in [NieOTj (see also [BNlOj and [Boull] for 
refinements), together with some other geometric arguments, should give 
the desired result. 
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